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Recent observations 1 2 3 have revealed an unexpectedly high 
binary fraction among the Trans-Neptunian Objects (TNOs) 
that populate the Kuiper Belt. The TNO binaries are strikingly 
different from asteroid binaries in four respects 2 : their frequency 
is an order of magnitude larger, the mass ratio of their compo- 
nents is closer to unity, and their orbits are wider and highly 
eccentric. Two explanations have been proposed for their for- 
mation, one assuming large numbers of massive bodies 4 , and one 
assuming large numbers of light bodies 5 . We argue that both as- 
sumptions are unwarranted, and we show how TNO binaries can 
be produced from a modest number of intermediate-mass bodies 
of the type predicted by the gravitational instability theory for 
the formation of planetesimals 6 . We start with a TNO binary 
population similar to the asteroid binary population, but sub- 
sequently modified by three-body exchange reactions, a process 
that is far more efficient in the Kuiper belt, because of the much 
smaller tidal perturbations by the Sun. Our mechanism can nat- 
urally account for all four characteristics that distinguish TNO 
binaries from main-belt asteroid binaries. 
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The TNO binary 1998WW31 has the following dynamical properties 2 : mass ratio 
1712/1711 ~ 0.7, eccentricity e ~ 0.8, semimajor axis a ~ 2 x 10 4 km, and inferred 
radii r\ ~ \.\ri ~ 10 2 km, where m\[r\) and 1712^2) are the masses (radii) of the 
primary and the secondary, hence ajr\ > 10 2 . All this stands in stark contrast to 
the properties of typical main belt asteroid binaries 7 , where m 2 /mi <1, e~0, and 
a/n < 10. 

Asteroid binaries are thought to have been formed by collisions 8 , in a scenario 
similar to the leading candidate for the formation of the Earth-Moon system 9 - 10 : 
two asteroids collide, leaving a small fraction of their combined matter with a large 
fraction of their relative angular momentum in a disk. Some of the disk matter then 
quickly coagulates into a small companion. The observed characteristics, m 2 /mi <C 1, 
e ~ 0, and ajr\ < 10, are all natural consequences of this scenario 11 . 

Clearly, we need a different mechanism for the formation of 1998WW31. We can 
look for an analogy by considering dynamical binary formation among stars in dense 
stellar systems, where there are three channels: 1) tidal capture, where tidal dissipa- 
tion during a close encounter between two single stars leaves the system bound 12 ; 2) 
three-body binary formation, in a simultaneous close encounter between three single 
stars where one of the three stars carries the excess energy away, leaving the other 
two stars bound 13 ; and 3) exchange reactions, where a single star encounters a binary, 
and replaces one of its original members 13 . 

Channel 1 is analogous to the standard scenario for asteroid binary formation. It 
will indeed occur: each TNO has grown through accretion, and much of this accretion 
has happened through collisions with an object comparable in mass to the growing 
TNO itself 14 ; in some of these collisions there will be too much angular momentum to 
form a single stable body, and in these cases it is unavoidable to form tight circular 
binaries of strongly unequal mass. Such binaries, when formed as intermediary stages 
during the growth of a TNO, are short-lived. The cross section for disruption of a 
binary is related to the cross section for accretion by the ratio of the sizes of the 
binary orbit and the TNO primary. Therefore, the binary will be broken up by a 
perturbation of a third body with mass comparable to the primary well before such 
a body will hit the primary itself, in the vast majority of cases. 

Channel 2 would require a near-simultaneous encounter of three massive objects 
with low enough velocities to allow an appreciable chance to leave two of the objects 
bound. For this to work, the random velocities of the most massive objects should 
be significantly lower than their Hill velocities. Under such conditions, this channel 
could play a role, as pointed out by Goldreich et al. 5 , who assumed that there are 
~ 10 5 100 km-sized object embedded in a sea of small (< 1km) objects. This as- 
sumption, however, is at odds with Goldreich and Ward's theory for the formation 
of planetesimals 6 through gravitational instability, and it is hard to see how objects 
in the Kuiper Belt could form from non-gravitational coagulation, because the time 
scales are far too long 15 . In contrast, the gravitational instability theory predicts the 
size of the initial bodies to be 10 — 100km. Starting with these larger bodies would 
make channel 2 ineffective, because the velocity dispersion would be higher than the 
Hill velocity 16 - 17 . 

Recently, Weidenschilling 4 proposed a variation on the idea of using interactions 
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between three unbound bodies in order to create a binary. He studied how a third 
massive body could capture the merger remnant from a collision of two massive 
bodies if the third body were near enough during the time of the collision. This 
mechanism seems unlikely to work, however, since it requires a number density of 
massive objects about two orders of magnitude higher than the value consistent with 
present observations 5 . 

Goldreich et al/ J have proposed another mechanism, based on the dynamical fric- 
tion from a sea of smaller bodies that can turn a hyperbolic encounter between two 
massive bodies into a bound orbit under favorable conditions. Effectively, this mech- 
anism makes use of a superposition of three-body encounters, since each light body 
interacts independently with the two heavier ones, and in that sense it is another 
variant on channel 2. As we mentioned above, the gravitational instability theory for 
the formation of planetesimals 6 would exclude the existence of such a sea of small 
objects, and since the alternative theory of nongravitational agglomeration does not 
seem to work, we will explore the consequences of dropping channel 2, which leaves 
us with channel 3. 

Channel 3 can only operate when there are already binaries available for encoun- 
ters. The only binaries that are expected to be formed frequently are the ones pro- 
duced in channel 1, so we should consider channel 1 and 3 in tandem. What remains 
to be done is to check whether the resulting properties of the binaries produced are 
the ones we are looking for in the Kuiper belt, and to check whether the formation 
efficiency is high enough to explain the abundant presence of surviving TNO binaries. 

Starting with the first task, consider a relatively massive TNO primary in a bi- 
nary orbit with a much less massive secondary, embedded in a sea of smaller particles, 
most of which are far less massive than either of the binary components. The small- 
est particles, when they come close, will simply accrete on either the secondary or 
primary or just pass through the system. The cumulative effect may drive the binary 
components to collide, which is fine, since then sooner or later a large impact is likely 
to create again a new temporary binary of a similar type. However, if the binary 
encounters a particle with a mass m that is comparable to the mass of the primary 
component (mi ~ m m 2 ), the most likely result is an exchange reaction, in which 
the incoming object replaces the original secondary 18 . Figure 1 shows an example of 
such a reaction. 

Using the impulse approximation for the replacement, together with conservation 
of energy and specific angular momentum, we can make a crude estimate for the 
semimajor axis a and eccentricity e of the newly formed binary in terms of the corre- 
sponding values ao and eo = for the original binary. Viewing the replacement as an 
isolated two-body encounter, we see immediately that the asymptotic escape speed of 
the initial secondary must be comparable to that of the incoming velocity at infinity of 
the third body. Since the latter is much smaller than the relative velocity of the binary 
components (as a necessary condition for runaway growth), the binding energy of the 
binary will not change much during the exchange, hence mim2/ao ~ m\mja where 
a is the new semimajor axis after the exchange. This implies a /a® ~ mjvrt2 3> 1: the 
size of the orbit of the binary is increased in proportion to the increase in mass of 
the object orbiting the primary. Under the impulse approximation, the interaction 
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Figure 1: An example of a binary-single-body exchange interaction, in the '(massive, 
light) meets massive' category discussed in this paper. Bodies 1 and 2 have masses 
mi — 1 and m 2 = 0.1, respectively, forming a binary with an initially circular orbit. 
Body 3, with mass = 1, encounters the binary on an initially parabolic orbit. In 
panel (a), the whole scattering process is shown. Panel (b) shows the complex central 
interaction in more detail, while panels (c) and (d) show the orbits of the initial and 
final binary, respectively. Note that the final binary orbit is highly eccentric and 
much wider than the initial circular binary orbit. 
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happens in a space small compared to the distance a to the primary. Neglecting 
factors of < 2 related to reduced mass, this implies that the angular momentum per 
unit mass of the newly bound particle with respect to the primary is the same as 
that of the previous companion. Conservation of specific angular momentum of the 
system then gives 1712^(1 — e ) ~ ma(l — e) which gives 1 — e m 2 /m <C 1: the 
eccentricity of the new binary is almost unity and the orbit is very elongated. 

While the impulse approximation may not be a very good assumption, at least 
qualitatively it predicts that exchange reactions provide us with just what we wanted: 
binaries with small mass ratios in large highly eccentric orbits. In order to perform 
a more quantitative check, we have run a series of scattering experiments to obtain 
the relevant cross sections. Here, we assumed that the initial binary has a mass ratio 
of 20 : 1 and semimajor axis a® = 20r\, where r\ is the radius of the primary. These 
values are typical for main-belt binary asteroids, with 1712/1711 < 0.1, and separations 
5 — 40 times the radius of the primary. We choose parabolic relative orbits for the 
single body approaching the binary, with periastron distances uniformly distributed 
between and 20a . In case of resonance reactions, where the incoming body is 
captured, and the whole system undergoes a complicated three-body dance, we only 
followed the system as long as all three bodies stayed within a maximum distance of 
lOOOao from the other two. If this condition was violated for any of the three bodies, 
we considered that body to escape, due to the perturbing tidal field of the Sun, which 
in the Kuiper belt corresponds to a Hill radius of order of 10 3 a . 

Table 1 shows the cross sections for those processes in which the initial binary 
membership is altered. Channels (a), (c) and (e) result in binaries with two mas- 
sive components, while channels (b) and (d) produce binaries with large mass ratios. 
There are four remaining channels that do not produce any binary: a triple merger 
1+2+3, and three channels in which two of the three bodies merge while the third 
body escapes. The cross sections of these four cases are summed together under 
(f). Note that about 80% of these interactions result in binaries with two massive 
components. These cross sections were calculated using a scattering code that incor- 
porated an effective tidal cut-off, and the results were checked independently through 
a comparison with the starlab three-body scattering package 19 . 

Figure 2 shows the normalized differential cross sections for the semimajor axis 
a and eccentricity e, for forming a final binary with two massive components. The 
distribution for the semi-major axis is strongly peaked at a = 20, in good agreement 
with the simple argument presented above. Similarly, the eccentricity peaks at 0.95, 
as expected. For the physical values of the semi-major axis, we have assumed a radius 
radius for the primary of r\ = 75km, the estimated size of the primary of 1998WW31, 
based on the mass deduced from the binary motion and an assumed mean density of 
lg/cm 3 . Figure 3 shows the distribution of the final binary in the a, e plane. The 
orbital elements of 1998WW31 are consistent with the binary having formed through 
the processes modeled here. 

We are now in a position to confront our second task: to check whether the ex- 
change channel is efficient enough to produce the observed binaries. A straightforward 
approach would be to derive semi-analytic estimates or to perform numerical simu- 
lations within the context of a model for the protoplanetary disk. Such an approach 
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Figure 2: Normalized differential cross sections for the formation of a 'massive- 
massive' binary, under the conditions specified in the text (channels a, c and e in 
table 1), with respect to the semi-major axis a (top panel), and eccentricity e (bottom 
panel) of the final binary. The initially circular binary has a = 1 in the dimensionless 
units used for da /da, while the physical units are given for reference at the top of the 
figure. The filled points are the total values for the differential cross sections, while 
the open circles are the contributions from the merger channels (c and e in table 1). 
Note the double-peaked structure in the top panel: the sharp peak toward a ~ 20 
arises from non-resonant exchanges, where the final binary has an energy comparable 
to that of the initial binary; the broad peak around a ~ 10 arises from resonant 
exchanges, where the memory of the initial binary is wiped out, leading on average 
to more strongly hyperbolic escape in which a harder binary is formed. 
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Figure 3: Orbital properties of 'massive-massive' binaries formed in our scattering 
experiments: a and e have the same meaning and units as in fig. 2. Contributions 
from exchange reactions, channel (a) in table 1, are limited by energy conservation to 
a < 20, and give rise to the horizontal rim in the middle of the figure. Contributions 
involving mergers, channels (c) and (e) in table 1, can lead to a values all the way 
to the Hill radius a 10 3 , but are limited by angular momentum conservation to 
increasingly high e for increasing a. The star symbol shows the observed orbit for 
1998WW31. Boxes around the star indicate the observational 1- and 2-a error bars. 
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Table 1: Cross sections a for various configuration-changing channels in binary-single- 
body scattering. The gravitational focusing factor v 2 is scaled out in order to obtain 
finite values in the parabolic limit, where v is the initial relative velocity between 
binary and single body at infinity. We use units in which G — m\ — m 3 = a = 1, 
where G is the gravitational constant, m\ and m 3 are the masses of the heaviest body 
in the binary and the single body, respectively, and a is the initial semi-major axis of 
the binary. The mass of the lighter body in the binary is m 2 = 0.05. The radii are 
r\ = r,3 = 0.05 and r 2 = ri^/mi) 1 / 3 m 0.01842. The scattering processes are coded 
as follows: (x, y) indicates a binary in the final state with components x and y, while 
a p + q indicates the product of a merger between bodies p and q. A single body 
z in the final state is indicated by (,),z. The physical meaning of the six channels 
is as follows: (a) an exchange reaction resulting in a massive-massive binary; (b) 
an exchange reaction resulting in a massive-light binary; (c) a merger resulting in a 
massive-massive binary; (d) a merger resulting in a twice-as-massive-light binary; (e) 
a merger resulting in a massive-massive binary; (f) no binary is left, after three-body 
merging or two-body merging followed by escape. 

was taken by Weidenschilling 4 within the context of simultaneous encounters of three 
massive bodies, and by Goldreich et al. 5 within the context of encounters of two mas- 
sive bodies while taking into account the dynamical friction from the background sea 
of lighter bodies with or without the additional passage of a third body. In both 
cases, they had to assume local disk densities high enough for their scenarios to work. 

Interestingly, for our scenario we can simply bypass such an approach, through 
the following observation. Starting with TNOs of intermediate mass, as predicted 
by Goldreich and Ward's theory for the formation of planetesimals 6 , the heaviest 
TNOs will accrete mass primarily through collisions with other TNOs of comparable 
mass 16 - 14 . We also know that many of these collisions are of the 'giant impact' type 
that give rise to the formation of a tight circular strongly unequal-mass binary as 
discussed above in channel 1. We therefore have to make no assumptions whatsoever 
about the density or other physical parameters of the protoplanetary disk, since 
changing these parameters will change the time scales for accretion and collisional 
binary formation in exactly the same way. All we have to do is estimate what fraction 
of encounters between comparable-mass TNOs will give rise to 'giant impact' type 
binaries, and how long such binaries survive on average before they are destroyed 
again. 

For example, let us assume that one in three collisions between comparable TNOs 
gives rise to a binary. This is not unreasonable, since in many collisions the encounter 
will be far from head-on, with an angular momentum that would exceed that of the 
merger product even if it would rotate at break-up speed, forcing the formation of 
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an accretion disk that is likely to form a small companion. In two of the three cases 
in which no binary is produced, we have to wait until another major collision occurs. 
Let us introduce this typical waiting time as T. What will happen, and how fast will 
it happen, in the remaining one out of three cases in which we form a binary? 

In the gravitational focusing regime, the cross section for significant three body 
interactions to occur is set by the size of the orbit a , while the cross section for 
direct collisions is set by the size r of the TNO primary. Therefore, our newly-formed 
binary is likely to undergo an exchange reaction on a time scale (r/ao)T, significantly 
smaller than T. As a result, a will increase significantly, as we have seen above. Strong 
three-body interactions will subsequently occur on an even much shorter time scale 
(r/a)T <C T. After the first exchange reaction, a very wide roughly equal mass binary 
exposed to encounters with a third body of comparable mass will on average shrink 
the orbit and randomize the eccentricity 13 . As a result, the semimajor axis will shrink 
systematically, while at any given time the eccentricity will still be considerable, since 
the 'thermal' distribution /(e) = 2e of eccentricities that follows phase space volume 
is tilted towards high eccentricities. 

When the orbit becomes small enough, subsequent three-body encounters are 
more and more likely to lead to a physical collision between two or three of the 
TNOs involved. If all three collide, we are back where we started, and the resulting 
system may be a single body (with an assumed chance of 2/3) or a strongly unequal- 
mass binary (chance 1/3). If two of the bodies collide, the third one may remain in 
orbit, with a thermal (e.g. high) eccentricity expectation, or it may escape. In the 
latter case, we again are back where we started. In the former case, we still have 
an equal-mass and likely highly eccentric binary. Even its separation will still be 
high, for the following reason. Long before the semimajor axis becomes comparable 
to the size of the bodies, a typical three-body resonance reaction is likely to lead 
to a collision between the three bodies, since each pass through the resonant period 
gives a renewed chance for a collision. The chance for collisions in resonant encounters 
becomes significant 20 when r/a ~ 0.03. For simplicity, let us assume that an exchange 
reaction turns a 'giant impact' binary into a binary with a semi-major axis of a ~ 
300r. Each subsequent strong encounter will on average decrease a by a factor 21 
~ 1.2. We thus need to wait for a dozen such encounters to occur before reaching 
a ~ 30r and facing a significant chance for a collision. The time scale for each 
encounter to occur is ~ (r/a)T. The waiting time for the last encounter in this series 
to occur, under these simplifying assumptions, is (1/30)T, while each previous waiting 
time was less by a factor 1.2. Summing this series, we get a total waiting time of 
(T/30)(l + (1/1.2) + (1/1.2) 2 + ...)« (T/30)/(l - (1/1.2)) = 0.2T, as an estimate 
for the duration for those dozen encounters to happen. 

Under these assumptions, in 1/3 of the cases, we wind up with an equal- mass 
TNO binary with the observed properties for a period ~ 0.2T, compared to a 2/3 
chance to wind up with a single TNO for a period ~ T. This allows us to derive the 
rate equation for the formation and destruction of the binaries. If we denote by Ns 
and Nb the number of single bodies and the number of binaries, respectively, we have 

^ = l Ns -± 2 -N B 
dt 3 0.2 3 
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dN s dN B 
~dT ~ ~ dt 

If we measure time in unit of T. So for the stationary state we have dN B /dt = 
dNs/dt = 0, and N B = 0.2N s /2 = 0.1N S . Therefore, the binary fraction is ~ 10%. 
When accretion in the Kuiper belt region diminished, the number of single and binary 
objects was frozen, with a ratio similar to this steady-state value. 

Clearly, this predicted binary fraction of 10% is dependent on the assumption we 
made for the size of a newly formed equal-mass binary, and our other estimates have 
also been rather crude. It is clear, however, that we expect to find any TNO as part 
of a very wide eccentric equal-mass binary during at least a few percent of its history. 
This implies that among all TNOs, after cessation of the accretion stage several 
percent or more were accidentally left in such a binary phase. The fact that more 
than 1% of the known TNOs are found to be in wide roughly equal-mass binaries 
is thus a natural consequence of any accretion model independent of the assumed 
parameters for the density and velocity dispersion of the protoplanetary disk or the 
duration of the accretion phase. 

We conclude that we have found a robust and in fact unavoidable way to produce 
the type of TNO binaries that have been found. As a corollary, we predict that 
future discoveries of TNO binaries will similarly show roughly equal masses, large 
separations, and high eccentricities. 
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